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One long-standing difficult problem in quantum dissipative dynamics is to solve the spin-boson 
model in a non-Markovian regime where a tractable systematic master equation does not exist. The 
spin-boson model is particularly important due to its crucial applications in quantum noise control 
and manipulation as well as its central role in developing quantum theories of open systems. Here 
we solve this important model by developing a non-Markovian quantum Langevin approach. By 
projecting the quantum Langevin equation onto the coherent states of the bath, we can derivie a 
set of non-Markovian quantum Bloch equations containing no explicit noise variables. This special 
feature offers a tremendous advantage over the existing stochastic Schrodinger equations in numerical 
simulations. The physical significance and generality of our approach are briefly discussed. 

PACS numbers: 03.65.Yz, 42.50.Lc, 05.40.Ca 


I. INTRODUCTION 

Quantum Langevin equation (QLE) provides a direct 
depiction of the temporal behaviors of physical observ¬ 
ables under the influence of a bath of quantum parti¬ 
cles As such, QLE has many important applica¬ 

tions in quantum optics 0, the input-output theory Q, 
and the quantum dynamics of dissipative atoms [8|- 
ini. For deriving a generic Langevin equation, however, 
Markovian approximation was usually employed to arrive 
at a tractable equation of motion. QLE beyond Marko¬ 
vian approximation can be also formulated to study the 
intriguing non-Markovian dynamics of dampe d q uantum 
systems and the Brownian motion systems IlM. A 
method [20l - [^ based on the Mori expansion mm can 
solve these kind of problems conveniently. The main idea 
of this method is to expand the time dependent operator 
with a set of time independent basis operators. This set 
of basis and the corresponding coefficients are govern by 
two recurrence relations. 

In the last decade, the so-called non-Markovian quan¬ 
tum state diffusion (QSD) equation [26l - [^ has been 
formulated nonperturbatively, so it can apply to the 
cases with strong coupl ing s between systems and envi¬ 
ronments (see, e.g., j29l - lMl |). The non-Markovian QSD 
has provided a powerful tool in numerically simulating 
many interesting physical models [13, 1131 • In particu¬ 
lar, high-order numerical methods for the non-Markovian 
QSD [3J-[3^ have been developed very recently, making 
some previously intractable problems becoming numeri¬ 
cally tractable. In fact, the QSD equation is a stochastic 
Schrodinger equation and it is solved by invoking the 
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noise realizations. For the important case of spin-boson 
model [131, however, matters are not as simple as the 
form of this model due to the fact that the spin-boson 
model does not admit an analytical treatment and an 
efficient numerical simulation is prohibited without in¬ 
cluding the higher-order perturbations [33 - l36j| . 

In this paper, we develop a stochastic quantum 
Langevin approach to solving non-Markovian quantum 
dynamics of the spin-boson model. This model which 
is the multi-mode case of the quantum Rabi model[23l- 
[Isjl involves non-conserving processes due to the counter¬ 
rotating terms. Consequently, it poses a long-standing 
difficult problem in studying non-Markovian quantum 
dynamics [33| . By projecting the non-Markovian QLE 
onto the coherent states of the bath, we convert the op¬ 
erator QLE into a c-number stochastic QLE, which is 
formally analogous to the non-Markovian QSD equation. 
Therefore, the useful techniques developed for the QSD 
can apply to the c-number stochastic QLE as well. Re¬ 
markably, we find that the stochastic QLE can be further 
reduced to a set of simple non-Markovian quantum Bloch 
equations without involving any noise variables. This 
provides a much more efficient method to solve the non- 
Markovian quantum dynamics of the spin-boson model. 
As shown below, the method developed here is quite gen¬ 
eral, so it may offer significant numerical advantages for 
simulating open quantum systems coupled to bosonic en¬ 
vironments when higher-order perturbation is unavoid¬ 
able. 

The paper is organized as follows. In Sec. II, we obtain 
a stochastic QLE by projecting the non-Markovian QLE 
onto the coherent states of the bosonic bath. Then, in 
Sec. Ill, we convert the stochastic QLE into a c-number 
stochastic QLE, which is formally analogous to the non- 
Markovian QSD equation. In Sec. IV, we further reduce 
the c-number stochastic QLE to a set of simple non- 
Markovian quantum Bloch equations. The Extensions 
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to the cases of complex correlation function and finite 
temperature are discussed in Sec. V and VI, respectively. 
Finally, Sec. VII gives the conclusion of our work. 


II. STOCHASTIC QLE 

The spin-boson model is described by Htot = 
with (setting h=l) 

k 

Hint = 'y ^ • (1) 

k 

Here Hq is the Hamiltonian of the uncoupled spin and 
multi-mode bosonic bath. Hint models the interaction 
between the spin and the bosonic bath, Gy and Gz 
are Pauli operators, and a\, {ak) is the fcth-mode bosonic 
creation (annihilation) operator of the bath. We assume 
that the state of the total system is initially factorized as 
l^'o) = |'0) ® |0), where the bosonic bath is in the vacuum 
state |0) (i.e., at zero temperature). 

The interaction Hamiltonian Hint can be rewritten as 
the sum of rotating and counter-rotating terms, 

Hint = ^(5fc4cr_ -H gkakCr+ + 5fc4cr+ + gkGkG-), 
k 

with Gx = G+ + G-. The counter-rotating terms g’^aj,G+ 
and gkGkG- break the conservation of excitation number, 
giving rise to high-order noise appearing in the stochastic 
equation of quantum dynamics |26l428l |. 

Starting from the Heisenberg equations of the Pauli 
operators and the bosonic operators of the bath, one ob¬ 
tains 


dt 

dt 

d 

dt 


Ckxit ) = - UJGy { t ), 

Gy(t) = ujGx{t) - 2 ^ Gz{t) gkak{t) + glalit) 


Gzit) = 2'^Gy{t) gkak{t) + glal{t) 


—ak{t) = -iukGkit) - iglGx{t). 
at 


( 2 ) 


where 


Gj{t) = j = X, y, z, 

are Pauli operators in the Heisenberg picture, and 
ak{t) = 

are the field operators in the Heisenberg picture. The 
field operators ak{t) and aj(,(t) in Eq. Q can be formally 
solved 


ak{t) = - igl f dse-“'=(*-*V,(s). 

( 3 ) 

Substituting the field operators in Eq. m with the formal 
solution in Eq. ([3]) we can achieve the QLE, 

^crx(t) = -U}Gy{t), 

^cry(t) = IXGxit) - 2Gz{t) [C(t) +^■'■(0] +Wz{t), 
^cr^(t) = 2Gy{t) [C(t) -f C^(t)] - Wy{t), (4) 

with a{t,s) = is the correlation func¬ 
tion of the bath, ^{t) = J2k defines a noise 

operator, and 

Wj{t) = 2iGj{t) [ ds[a{t,s) - a*{t,s)]Gx{s), 

Jo 

with j = y,z. In this QLE, both ^{t) and ^^t) act as 
“random noises” acting on the spin. 

Here we first consider the case of real correlation func¬ 
tion a{t,s) = a*(t,s), so that Wy{t) = Wz{t) = 0 in 
Eq. (HI). Note that the typical Ornstein-Uhlenbeck cor¬ 
relation function a{t,s) = is indeed a real 

function. We define Bargmann coherent states for the 
bosonic bath. 


z) = (^ \zk) = (5) 

k 


which satisfy 


ak\z} = Zk\z), al\z) 



( 6 ) 


When projected onto the Bargmann coherent states, the 
QLE in Eq. (|4|) is then converted to 


d . ^ 


d . ^ 


-ujGy{t;z), 

ujGx{t;z) 


-2 


+ 


2 



dsa[t, s )-— 

0 Zs 


dsait, s )-— 
bZs 


(^z{t]z), 

Gv{t]z). 


( 7 ) 


This is a stochastic QLE with the noise Zt = 
Sfc5fee““'=*2;fc- In Eq. (j^, Gj{t;z) = Gj{t)\z){z\, with 
j = X, y, z, and the functional chain rule. 


d _ f , dzs S 

dzk ~ J dzkSz/ 
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is used. Note that 

= n = 7W {crj(t;2)} . (8) 

fc ^ 

When statistically averaging Eq. m over all noise vari¬ 
ables via Eq. ([5]), one can recover Eq. m back to the 
QLE in Eq. (U). 


III. C-NUMBER STOCHASTIC QLE 


To convert the stochastic equation of operators [i.e., 
Eq. 0] into a c-number equation, we introduce the ex¬ 
pectation value of an operator a as (tr) = (4'o|cr|'I>o). 
Then, we have 


— {a^{t]z)) = -uj{ay{t]z)), 

d 

^{<Jy{t;z)) = uj{(T^{t;z)) 


_9 

m 


{a^{t;z)) = 2 


zt+ [ dsa(t,s)-^ 

Jo 

zt-\- dsa{t,s) 

Jo 


{^zit;z)), 


_S_ 

Szs 


(cry(t;z)). (9) 


we can write the c-number stochastic QLE in a time-local 
form, 


dt 

where 


A{t, z) = [-iV. + Czt + CO{t, z)] A{t, z), (13) 


0{t,z)= / dsa{t,s)0{t,s,z). 


Also, using Eq. dT51) and the relation 

S d . d S . 


we obtain the equation for 0{t,s, z) operator. 


_9 

m 


0{t, s, z 


) = [-ili + Czt -I- CO{t, z), 0{t, s, z)] 


60{t,z) 

SZs 


(14) 


As in Ref. [2^ the initial condition of the 0{t, s, z) oper¬ 
ator can be derived as 0{t,t,z) = £. 


IV. NON-MARKOVIAN QUANTUM BLOCH 
EQUATION 


Define A{t,z) = z)), {ay{t-, z)), {az{t; z))]'^, where 

T denotes the transpose of a matrix. Equation m can 
be written in a matrix form as 

d 

— z) = z) 4- CztA{t, z) 

ot 

+C / dsa{t,s)- — A{t,z), (10) 

Jo ^Zs 

with 


(^ 

—ioj 

0\ 


/o 

0 


iuj 

0 

0 

, C = 

0 

0 

-2 

V 0 

0 

0/ 


lo 

2 

oj 


Formally, this c-number stochastic QLE is analogous 
to the non-Markovian QSD equation governing non- 
Markovian quantum trajectories . The difference here 
is that the QSD equation is a stochastic differential equa¬ 
tion for quantum states of the system (i.e., a stochas¬ 
tic Schrodinger equation), while the c-number stochastic 
QLE in Eq. (fT(Tl) corresponds to a stochastic differential 
equation of physical variables. 

Here we introduce 0{t, s, z) operator by 

z) = 0{t, s, z)A{t, z). (12) 

uZg 

Note that although we use the notation 0(t,s,z) which 
is similar to the O operator used in QSD approach, their 
meanings are different. Here 0{t,s,z) is defined for an 
arbitrary operator, rather than for a quantum state. Now 


To obtain the desired quantity 
A{t) = M{A{t, z)} 

(15) 

where {(Jj{t)) = {Ai{aj{t; z)}), one can numerically solve 
Eq. (1131) for each realization of the noise Zt and then im¬ 
plement the statistical average, as in the case of numeri¬ 
cally solving QSD equation. However, when higher-order 
perturbation is involved in the QSD, one must pay the 
price of long computation time in order to achieve accu¬ 
rate results. Below we show that with our QLE approach, 
this simulation process can be significantly sped up. 

By directly implementing statistical average on 
Eq. (fT^ . we have 

^A{t) = —iT-LA^t) + CJA {ztA{t, z)} 
ot 

+CM{d(t,z)A{t,z)}. (16) 

In Eq. (fTBl) . Ai {ztA{t,z)} can be written as 

M{ztAit,z)} = llJ ^e-l^'=l4t(vI/o|H(t)k)(^|'I/o), 

where B{t) = {aj;{t),(Ty{t),az{t))'^. Because zt\z) = 
^{t)\z), we have 

M{ztA{t,z)} = {^o\Bm{t) 

k 

= i'^oiBimmo) 

= 0, (17) 
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where we have used the relation ^(t)|4'o) = 

| 0 ) = 0 . 

It is known that the 0{t,s,z) operator can be ex¬ 
panded as 


0{t,s,z) = Oo(t, s) -k 



X Zv^ . . . Zy^dvi . . .dVn- (18) 


Because 


M {dn{t,Vi,...,Vn)Zy^ . . . 2:„„(4'o|S(i)|-z)(-2|^o)} 

= M {dnit,vi,... ,Vn){'^o\B{t)({vi )... |z)(z|4'o)} 

= On{t,Vi, . . . ,Vn){'^o\B{t)^{vi) . . . ^(Vn)M {|z)(z|}4'o) 
= On(t,Vi, . . . ,V„)('i>olB(t)((vi) . . .C(Wn)l^'o) 

= 0, (19) 

where 


On(t,Vi, ... ,v„) = / dsa(t,s)On(t,s,vi,...,v„), 

do 

Eq. (fTO) is finally reduced to our central result 

^A(t) = -mA(t) -k £do(t)A(t). (20) 

ot 

Here we call it a non-Markovian quantum Bloch equa¬ 
tion, in which no noise variables are involved. Now note 
that only the noiseless term of the functional expansion 
in Eq. (jl 8 p is important in solving the non-Markovian 
quantum dynamics of the system. Because no noise vari¬ 
ables are involved, Eq. can be numerically solved 
very efficiently. 

Figure 1 shows the time evolution of (cr^) for a bath 
with the Ornstein-Uhlenbeck correlation function. It can 
be seen from Fig. 1(a) that when increasing 7 , the envi¬ 
ronmental memory time I /7 decreases and (cz) exhibits 
a clear transition from an oscillation to an exponential 
decay. Physically, this is to some extent connected to 
the overdamped oscillator, where increasing the friction 
on the velocity of the oscillator has the effect of turning 
an oscillation into an exponential decay. 

The results are obtained by solving the non-Markovian 
quantum Bloch equation (l20l) . with Oo{t) determined by 
(see Appendix A) 

^Ooit) = -i[n,do{t)\ + [cdo{t),do{t)\--ido{t) 

Ey 

+^£ + /:Qi, 

d ^ 

Q„] + ^[£Qfc, Qn-k] - in + 1)7Q„ 

+ -^[£j Qn-l] + (’T'+ l)£Qn-|-l, (21) 

with initial condition Qq = Oo{t). These hierarchical 
equations do not contain any explicit noise variables. In 



a)t 

FIG. 1. (color online) Time evolution of {aA for a bath with 
Ornstein-Uhlenbeck correlation function a(t, s) = 

(a) Af = 100 , and the inverse of the correlation time is chosen 
to be 7 = 0.2, 0.4 and 0.8, respectively, (b) 7 = 0.2, and 
the hierarchical order is chosen to he Af — 0, 3, 10 and 100, 
respectively. Also, the coupling strength is chosen to be yF = 
0.2 in both (a) and (b). 


numerical calculations, one can truncate Eq. (EB) at a 
given hierarchical order M by choosing Q^f+i = 0. The 
results in Fig. 1(a) are very similar to those in Ref. [s^ 
obtained using the QSD method, showing apparent non- 
Markovian behaviors at small values of 7 . In [s^, the 
simulations for the curve with 7 = 0.2 took about 36 
days to execute on an Intel core-i7 CPU core, but only 
a few seconds here by solving the non-Markovian quan¬ 
tum Bloch equation (l20l) via the noiseless hierarchical 
equations in Eq. (EB). This is because of the numeri¬ 
cal efficiency of our method without invoking any noise 
realizations. While {aA at A/” = 0 and 3 deviate from 
those at A/” = 10 and 100 , (ctz) at A/” = 10 and 100 look 
nearly identical [see Fig. 1(b)], revealing fast convergence 
of our results with the hierarchical order Af. In contrast, 
the results of (cr^) obtained using the QSD method have 
considerable differences between the A/” = 10 and 100 
orders of the hierarchical equation (see Fig. 2 in M), 
indicating much slower convergence with Af there. 

Many years ago, a proposalwasmade to convert aQLE 
with correlated fluctuations into a set of coupled equa¬ 
tions [2fll - [2^ . It was originally developed to study the 
QLE with an additive noise (e.g., the quantum Brow¬ 
nian motion is such a case) and then extended to the 
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multiplicative-noise case [2lj. Here we study quantum 
dynamics of the spin-boson model. This model involves 
a multiplicative noise in the QLE and is a more com¬ 
plex, open problem of quantum statistical physics. The 
central point of our approach is to reduce the QLE to a 
simple differential equation with no noise variables, i.e., 
the quantum Bloch equation in Eq. (20). Moreover, we 
obtain a set of coupled equations, as in Refs. [20l - l^ . 
and then use it to efficiently calculate C>o(t) in Eq. (20) 
without invoking any noise realizations. This is the key 
reason our approach has a high numerical efficiency. 

V. EXTENSION TO THE CASE OF COMPLEX 
CORRELATION FUNCTION 

When the correlation function is complex, i.e., 
a{t,s) ^ a*{t,z), an extra term yV{t,z) = 

(0, {wz{t)\z){z\), {wy{t)\z){z\))’^ is added to Eq. ([131): 


temperature bath, where a fictitious bath with Hamil¬ 
tonian Hb = '^k{—^k)b\bk is introduced. The corre¬ 
sponding Hamiltonian of the total system then reads 

H = ^o-2+X! + 9ka}^ Wfc (^alak - b\b,^ . 

k k 

(26) 

When applying a Bogoliubov transformation [3^ to the 
system, 

ak = \/nk + Icfc -I- y/n^d\, 

bk — ^^k T Id/j; TlkC^ , (27) 

where hfc = — 1 ]“^, the composite bath of 

bosonic operators and 6 ^ initially prepared in a ther¬ 
mal state is equivalently converted to a virtual composite 
bath of bosonic operators Cfc and dk in the vacuum state 
|0) = |0)c ® \0)d, with Cfe|0)c = 0 and dk\0)d = 0. Now, 
the Hamiltonian of the total system is transformed to 


^M(t, z) = [-iH + Czt + CO{t, z)] A{t, z) + z). 

( 22 ) 

As the simplest approximation, one can apply a Marko¬ 
vian approximation only to the term W(t, z) in Eq. (1221) 
by taking a{t,s) in Wj{t) as a{t,s) = S(t — s). Then, 
W{t, z) = 0, and both the same c-number stochastic QLE 
(11311 and the same non-Markovian quantum Bloch equa¬ 
tion ( 1 ^ are thus obtained. 

Also, we can replace C7x(s) in Wj(t) by (Jx(t). Then, 

Wz(t) ~ —iv(t)ay{t), and Wy(t) ~ iv(t)a^(t), with v(t) = 

4 Jp dsIm{Q;(t, s)}. This approximation can give very ac¬ 
curate results at the early stage of quantum evolution. 
The term yV{t, z) in Eq. (1221) can be written as 

Wit,z) = -iV{t)Ait,z), (23) 

with 

/O 0 0 \ 

V(t) = 0 v{t) 0 . (24) 

V 0 0 -v{t) J 

Thus, Eq. (1^ is reduced to 

^M(t, z) = + Czt + CO{t, z)] A{t, z), (25) 

which has the same form as Eq. m, with only TL re¬ 
placed by 7^(t) =71+ V{t). Also, we can derive the 
equation for 0{t, s, z) operator and the non-Markovian 
quantum Bloch equation, which have the same forms as 
Eqs. (fT3)l and (l 20 |) . respectively, but with H replaced by 
'H(t) =71 + V{t) as well. 


VI. FINITE-TEMPERATURE EXTENSION 


H = ^ Vnk + IcTx -I- gkCk^ + ^ Wfcc[,Cfc 

k k 

+ ^ ^ yj^k^x 9k 4)-E uJkdldk- (28) 

k k 

Similar to Eq. o, the Pauli operators obeys the QLE 
^CTx(t) = -WCry(t), 

^cry(t) = - 2(j^{t) irit) + Ct(^) + WTz{t), 

^az{t) =2ay{t) ^t(<)+Ct(0 -WTy{t), (29) 

with the temperature-dependent noise operator 
?T(t) = ^ [Vnk + lgfce““'“‘cfc(0) -k V^5fce*‘^'“‘4(0)] , 


and 

= 2iaj{t) / ds[aT{t,s) - aT{t,s)]ax{s), 

Jo 

where j = y,z, and 

arit, s)=J2 \9k\" \i^k + 


is the finite-temperature bath correlation function. 

Using a similar procedures above, we can derive the 
c-number stochastic QLE at a finite-temperature as 




-4(t, x) = [-iJb -k Cxt + CO(t, x)] A(t, x) + WT(t, x), 

(30) 

where 


With the thermo-field method [H HI, [s^, we can Xt = ^ [s/nk + l 5 fce"*‘^'“*Zfc -k y/^glA^’^^Wk] 

map the finite-temperature bath onto a larger zero- 
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is the temperature-dependent noise, and 

Oit,x)=[ dsaT{t,s)0{t,s,x)- 
Jo 

The term Wt in Eq. (|3Q|) is 

yVT{t,x) = {0,{wTz{t)\zw){zw\),{wTy{t)\zw){zw\))'^, 

where \zw) = | 2 :){^|w), with the Bargmann coherent 
states defined by 

k) = l-Z/c) = 6 ^*= | 0 )c, 

k 

|w) = (^ |wfc) = (31) 

k 

which satisfy Ck\z) = Zk\z) and dk\w) = Wk\w), respec¬ 
tively. Note that Eq. (1501) is formally similar to Eqs. m 
and (|22l) . Therefore, we can solve the finite-temperature 
problem in an analogous way. 


VII. CONCLUSION 

We have developed a quantum Langevin approach to 
solving non-Markovian quantum dynamics of the spin- 
boson model. Instead of directly attacking the spin-boson 


model with our non-Markovian QLE, we arrive at a c- 
number stochastic QLE through projecting the operator 
QLE onto the coherent states of the bath. Furthermore, 
we have shown that the stochastic QLE can be reduced 
to a non-Markovian quantum Bloch equation. With the 
noiseless quantum Bloch equation, we can efficiently solve 
the non-Markovian quantum dynamics of the spin-boson 
model. In addition, we show that our approach is gen¬ 
eral enough to include the finite-temperature bath. Since 
the spin-boson model does not admit a non-Markovian 
master equation, therefore, generally one cannot arrive 
at a set of useful Bloch equations desirable from our ex¬ 
perience in dealing with Markov systems. We show in 
this paper that QLE paves a new avenue to bypass the 
stringent difficulty in deriving the non-Markovian master 
equations. We expect our stochastic quantum Langevin 
approach can play an important role for many other open 
quantum systems. 
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Appendix A: Derivation of the hierarchical equations in Eq. (|21|) 


From the equation of 0{t, s, z) operator in Eq. (fT4l) . it was obtained that the On(t, Vn) operators in 

Eq. (flSll obey the following hierarchical equation: 


Wt 


On{t,S,Vi, ...,Vn) = . .,?;„)] + [n + 1) COn+l{t, S, Vi, . . . , Vn) 


4-f E '^i-^^k{t,Vp^^i),...,Vp^^k)),On-k{t,S,Vp^^k+l),---,Vp^{n)], (AI) 

p,^es„ k=o 

with Oo{t,t) = C, On{t,t,vi, = 0, and On{t,s,t,vi ,... ,w„_i) = ^[C,On-i{t, s,Vi ,... ,u„_i)] for n > I. Here 

Sn denotes the permutation of all P„(fc)’s and (5„(t, vi,..., u„) = dsa(t, s)On(t, s,vi,..., Vn)- 
Let us define an operator 

rt 


[ ds 

f dv ,.. 

■ f 

lo 

lo 

Jo 


Qn{t) = / ds dvi... dVnOn{t,S,Vi,. . . ,Vn)a{t,s)a{t,Vi) . . .a{t,Vn), 

Jo Jo Jo 


(A2) 


with Qoit) = Oo{t)^ and consider a noise characterized by the Ornstein-Uhlenbeck correlation function a{t,s) = 
£2g-7 |t-s|^ It can be derived that 


dt 


Oo{t) =-^ J dsOo{t,s)a{t,s) 

a{t, s) + Oo(t,t)a(t,t) — X f dsOo{t, s)a(t^ s) 

Jo 


= ds 


= ds 






a{t,s) + —L-xOQ{t). 


(A3) 
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From Eq. dsa, it follows that ■^Oo{t,s) = —[i'H,Oo(t,s)] + COi{t,s) + [COo(t),Oo{t, s)]. Substituting it into 
Eq. (EH), we have 


^ao(t) = -[m,do{t)] + [c6o{t),do{t)]-^do{t) + ^c + CQi{t). 


This is the first equation in Eq. (EH). 
For n > 1, it can be derived that 


dt dt 


= ds dvi ... dV' 
Jo Jo Jo 

rt 


[ ds 

fdvi.. 

./ 

lo 

fo 

Jo 


■ 5 ^n) 


a{t,s)a{t,vi).. .a{t,Vn) 


+a{t,t) dvi... dVn-lOn{t,t,Vi,. .. ,Vn-l)a{t,Vi) .. .a{t,Vn-l) 

Jo Jo 

+na{t,t) ds dvi... dvn-iOn{t, s,t,vi,... ,Vn-i)a{t, s)a{t,vi)... a{t,Vn-i) 
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(A4) 


(A5) 


where we have used the relations On{t, t,vi,..., n„) = 0, and On{t, s, t,vi,..., n„_i) = On-i{t, s,vi,..., n„_i)]. 
Substituting s,vi,. .., Vn) in Eq. (lAll) into Eq. (IA5I) . we then obtain 


d ^ r 

Qn = —ilTd, Qn] + Qn-k] — {u + l)7Qn + [£, Qn-l] + (?^ + l)/iQn+l, (A.6) 

k=0 


dt 


which is the second equation in Eq. (EH- 
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